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In this paper we review the experimental and observational searches for stable, 
massive, elementary particles other than the electron and proton. The particles may 
be neutral, may have unit charge or may have fractional charge. They may interact 
through the strong, electromagnetic, weak or gravitational forces or through some 
unknown force. The purpose of this review is to provide a guide for future searches 
- what is known, what is not known, and what appear to be the most fruitful areas 
for new searches. A variety of experimental and observational methods such as 
accelerator experiments, cosmic ray studies, searches for exotic particles in bulk 
matter and searches using astrophysical observations is included in this review. 
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1 Introduction and General Considerations 



1.1 Scope of the review 

In this paper we review tlie experimental and observational limits on the existence of 
stable, massive, elementary particles - particles at least as massive as the electron, 
0.5 MeV/c^. We use the term elementary in a broad sense so that the presently 
known stable, massive, elementary particles are the electron - a basic elementary 
particle - and the proton - a composite elementary particle. 

The particles may be neutral, may have unit charge or may have fractional 
charge. They may interact through the strong, electromagnetic, or weak force, or 
through some unknown force, as long as the force allows direct or indirect detection 
of the particle. This review does not include the particle classes of neutrinos, axions 
or monopoles. 

As the reader already knows, except for the electron and proton, no stable, 
massive particles that fit our criteria have been found, although many have been 
proposed and sought. Our purpose is to provide a guide for future searches - what 
is known, what is not known, and what appear to be the most fruitful areas for new 
searches. Therefore this review includes a variety of experimental and observational 
methods: accelerator experiments, cosmic ray studies, searches for exotic particles in 
bulk matter, and astrophysical searches for signals from the annihilation of particle- 
antiparticle pairs. 

We have chosen the electron mass as the lower mass limit for this review for two 
reasons. First, in general, search methods differ above and below about 1 MeV/c^. 
Second, our own experimental interests are in more massive particles, in contrast to 
neutrinos and axions. 

Particle stability is a loose criterion for us. For example, we are interested in 
particles whose lifetimes are of the order of the lifetime of the universe so that 
searches for particles produced in the early universe can be meaningful. However, 
we are also interested in searches at accelerators for new particles with lifetimes 
sufficiently long to be directly detected, that is, longer than about 10~^ s. Finally, 
we do not classify bound quarks as stable particles, but we would classify as stable 
particles free quarks that lived long enough for direct detection. Thus in this review, 
a stable particle is operationally defined as one that can be observed in isolation 
and has sufficient lifetime to be detected by the search techniques described in this 
paper. 

While nuclei are not elementary particles, stable nuclei fit our operational search 
criteria. In this review we have not explicitly discussed searches for unknown, mas- 
sive, stable nuclei, however some of the search limits can be applied to the question 
of the existence of unknown, massive, stable nuclei. 
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We have kept this review to a reasonable length by imposing two drastic limita- 
tions. First, we give very few details of experimental and observational techniques. 
Second, while we note and sometimes summarize relevant particle and astrophysical 
theory, we do not give derivations. In writing a review paper the question always 
occurs as to the extent of the references. Here also we have chosen to be economical 
of space. References arc given for all quoted searches and for some reviews but 
we do not give references to general concepts in particle physics or astrophysics, 
concepts such as supersymmetry theory, early universe cosmology, or cosmic ray 
phenomenology. 

1.2 Dark matter 

Many present searches for massive elementary particles are motivated by the desire 
to identify the elementary particle or particles that compose what we call dark 
matter. Dark matter is thought not to be made up of baryons or charged leptons, 
possibly to be made up in part of neutrinos, but mostly to be made up of one or 
more types of unknown elementary particles. This review includes the results of 
searches for the major classes of proposed dark matter massive candidates, such as 
WIMPS (weakly interacting massive particles), SIMPS (strongly interacting massive 
particles) and CHAMPS (charged massive particles). 

Our search interests are more general. An unknown, stable, massive particle 
may exist in numbers much too small to explain the amount of dark matter in 
the universe. Still, the existence of such a particle would be of great importance. 
Indeed, one purpose of this review is to urge that a massive particle search effort 
not cease when the upper limit on that particle's density falls below dark matter 
density requirements. 

1.3 Massive particle production in the early universe 

In models for massive particle production in the early universe, models often moti- 
vated by the search for a solution to the dark matter puzzle, it is usually assumed 
the massive particle X and its antiparticle X, existed in thermodynamic equilibrium 
in the very early universe. Then, as the universe cools and expands, the X particles 
'freeze-out' of equilibrium and become thermal relics, their density depending upon 
their mass, Mx and their annihilation cross section (Jann- 

Since the X particles are non-relativistic, most of aann is in the S wave, leading 
to aann inversely proportional to M|-. The 'freeze-out' density, px is inversely pro- 
portional to aanm heuce proportional to M|-. Therefore the larger Mx, the larger 
Px, but Px cannot exceed the critical mass density of the universe. This leads to an 
upper limit on the mass of thermal relics so produced. Griest and Kamionkowski [1] 
give the upper limit of 3.4 x 10^ GeV/c^, for convenience we use 
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■mal— relic 

< 10^ GeV/c2. 



(1) 



There are proposals for getting around this upper hmit on Mthermai-reiic- For 
example Kolb et al. [2] have proposed processes by which particles with masses in 
the range of 10^^ to 10^^ GeV/c^ might be produced in sufficient quantity to explain 
the abundance of dark matter. If one is searching for massive particles with less 
than dark matter abundance, these non-thermal relic processes become even more 
attractive. 

Experimenters who look for massive particles with Mx > 10® GeV/c^ produced 
in the early universe must balance the time, work, and money necessary for such 
searches, against how much they beheve the proposals for getting around Eq. 1. 

1.4 Galaxy formation and halo particles 

During the formation of our galaxy there was an early period when the gravitational 
force of the galaxy mass controlled the velocity distribution of all individual particles. 
All particles have the same average velocity, Vavei calculated by using the virial 
theorem, namely: 

Vave ~ 300 km/s or (3ave ~ 10~^ (2) 

Here (3 — v/c. 

As the galaxy cools, the particle's kinetic energy decreases through coUisions 
with other particles. For an electron, proton or nucleus the kinetic energy is small 
and most electrons, protons and nuclei condense into ordinary matter. But there 
are two conditions that must be satisfied so that massive X particles also condense 
into, or are trapped by, ordinary matter: 

a. The interaction of the X with ordinary matter must be strong enough so 
that the X particle loses its initial high velocity through its collisions with ordinary 
matter. Electromagnetic interactions or perhaps strong interactions are required. 
An X particle with only the weak interaction may not be able to lose enough energy 
through coUisions. An example is the WIMP model where WIMPS with masses in 
the range of 50 to 500 GeV/c^ are assumed to still retain Pave ~ 10~^. 

b. If Mx is very large, say 10^^ GeV/c^, then its initial energy is 10® GcV and 
the particle will very rarely be trapped in ordinary matter unless it directly hits a 
large object such as the earth. In that case, depending upon the type of interaction, 
the X may be stopped in the earth's atmosphere or in the earth's crust or deep in 
the earth. Of course with sufficient energy the X might pass right through the earth. 

If massive particles of some type are not trapped by ordinary matter and remain 
moving in random orbits in the galaxy, we call them halo particles. 
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The mass boundary between when particles are trapped in ordinary matter and 
when they remain in the halo of the galaxy depends upon the interactions of the 
particle. The calculation of the mass boundary is not precise because of the different 
models used for particle interactions as the galaxy cools. Thus for a positive, unit- 
charge, massive particle, De Rujula et al. [3] give the boundary as 2 x 10^ GeV/c^, 
that is, above 2 x 10'^ GeV/c^ these massive particles remain in the halo. On the 
other hand, Dimopoulos et al. [4] give a much higher mass boundary, 10^ GeV/c^. 

1.5 Massive particle searches using colliders and accelera- 
tors 

Searches for stable, massive particles at high energy colliders - e+ e^, e''^ p, and 
p+ p^ - have two obvious advantages. First, the search mass range is known and, if 
a production cross section is assumed, the existence of a hypothetical particle may 
be directly tested. A simple example is the search for a stable, charged lepton, L='=, 
using 

e+ + e" ^ L+ + (3) 

A second advantage is that the stability requirement is less restrictive compared 
to other search methods, usually a hfetime greater than 10~^ s is sufficient for direct 
detection of the particle. 

The limitations of collider searches are obvious. The mass range is limited by 
the available energy. The upper limit to the Tevatron mass search range is several 
hundred GeV/c^ and that for LEP2 is about 100 GeV/c^. The significance of a null 
coUider or accelerator search is limited by our knowledge of the production cross 
section. Thus a search for free quarks cannot use conventional quantum chromody- 
namics (with its confinement hypothesis) to calculate the expected production cross 
section. 

There have also been some interesting studies using a fixed target accelerator. 
For example Prinz et al. [5] have searched for particles with very small charges. 

1.6 Searches in cosmic rays and halo particles 

There are three classes of particles impinging on the earth. The two observationally 
dominant classes are the primary cosmic rays coming from outside the earth, mostly 
proton and nucleons, and the secondary particles produced by interactions in the 
atmosphere such as muons and pions. The primary cosmic ray class is of more 
interest particularly because of the detection of very high energy particles, those 
with energies greater than 10^^ eV. Historically the class of secondary particles has 
been of great importance because this class was the source of the discoveries of 
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many particles. But. with the steadily increasing energy and intensity of colliders 
the secondary class is of decreasing interest. 

The third class of particles impinging on the earth consists of the assumed halo 
particles of the galaxy. Interest in this class has, of course, been greatly stimulated 
by the search for the constituents of dark matter. 

1.7 Searches for stable, massive particles in ordinary matter 

Stable, massive particles might be introduced into ordinary matter through two 
different processes. In one process occurring in the course of the cooling of the 
galaxy, massive particles might condense with baryonic matter to form stars, planets 
and smaller bodies such as asteroids. As discussed in Sec. 1.4, this process will occur 
if the particles are not too massive or if in their interaction with ordinary matter 
they have exceptionally large dE/dx. 

The other process that could introduce massive particles into ordinary matter 
would take place once stars and planets are formed. Some types of massive cosmic 
ray particles or halo particles impinging on the these bodies could lose sufficient 
energy to become trapped in the ordinary matter of these bodies. An example is 
the hypothesis that a massive, positive charge particle, X"*" falling through water 
might come to rest and form the exotic heavy water molecule HXO. 

2 Searches for Integer Charge, Stable, Massive 
Particles 

2.1 General considerations 

In this section we consider bounds on the existence of integer charge massive parti- 
cles, X. Wc assume that the particles can be positively or negatively charged and, 
although in principle the magnitude could be any integer, we restrict ourselves here 
to the case of unit charge. 

Limits on the existence of these particles come from a variety of experiments. 
Accelerator searches (Sec. 2.2) have a hmited mass reach but are able to detect 
particles with very small lifetimes. Experiments could also detect present day relic 
abundance of the X particles from production in some early epoch of the universe. 
Their subsequent history, which would tell us how best to search for them, depends 
on their mass and on theories of galaxy formation, as discussed in Sec. 1.4 (see also 
Refs. 3,4). The first of two extreme possibilities is that they remain in the galactic 
halo maintaining dynamical properties similar to those associated with cold dark 
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matter. The second possibility is that they condense out into the galactic disk and 
subsequently become trapped within ordinary matter. Searches in cosmic radiation 
probe the former scenario (Sec. 2.3) and searches in bulk matter (Sec. 2.4) test 
the latter. In addition, cosmic ray searches can also look for X particles produced 
and accelerated by high energy astrophysical processes. Finally, indirect limits on 
X particle properties can be inferred from a number of astronomical observations, 
as discussed in Sec. 2.5. 

Although there has always been interest in searches for trace amounts of particles 
with characteristics like that of X, the possibility that charged particles could have 
an abundance large enough to make them viable candidates for dark matter would 
at first thoughts appear ludicrous. However, as first shown by Dc Riijula et al. [3] 
and Dimopoulos et al. [4], if the particles are sufficiently massive they can indeed 
play the role of cold dark matter. This realization led to a revival of interest in 
searches for charged massive stable particles and much of the work reported below 
was a consequence of Refs. 3,4. To set the stage for later discussion (Sees. 2.3-2.5), 
we list here a few properties of dark matter that arc generally accepted (sec Primack 
et al. [6] for additional details). The first is the estimate of the dark matter density in 
the local solar neighborhood, pDM,iocai ~ 0.3 GeV/cm^. The second is the velocity 
distribution of dark matter, which is approximately Maxwellian with a cutoff at 
about 640 km/s (the local escape velocity) and with an rms velocity of about 260 
km/s. Thus, the flux of dark matter at earth can be approximated by: 

$ f« (0.3 GeV/cm^)(10-3c)/MoM ~ lO'^cm-h''^ /{MoM/GeV /c'^). (4) 

Here Mdm is the mass of a dark matter particle. Note that this flux is larger than 
the cosmic ray proton flux which is about [7] 1 cm~^s~^ for Mdm—^O^ GeV/c^. 
Finally, an upper limit on the density of dark matter is believed to be about 30% 
of the critical density of the universe [8], and about ten times more than that of 
baryons. On scales of our Galaxy, the dark matter resides in a much larger volume 
(the halo) than the baryons (which are constrained to the disk), hence the low value 
of 0.3 GeV/cm^ given above for the local neighborhood. 

Finally, as will be shown below, most experimental limits indicate that charged 
massive particles do not fill the budget demands for the density of dark matter. This 
has led to diminished interest in further searches for charged stable massive particles. 
However, as recent results are making clear, the energy budget of the universe is not 
dominated by one species, but is an amalgamation of baryons, neutrinos, exotic dark 
matter, exotic dark energy, etc [9]. Could it not also be that the dark matter itself is 
composed of a spectrum of different particles? We leave it to future experimenters 
to devote their efforts to searching for X particles whose abundance may be very 
rare in comparison to known species. 
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2.2 Accelerator searches for integer charge, stable, massive 
particles 

Equation 3, the pair production of integer charge particles in electron-positron an- 
nihilation, provides the most straightforward way to search for such stable massive 
particles. The beam energy and the measured momentum provide a determination 
of the mass of the sought particle, the pair production cross section is known and the 
LEP2 electron positron collider had a mass search range up to about 100 GeV/c^. 
Figure 1 shows the results of a search by Acciarri et al. [10], the L3 experiment 
at LEP. Up to the search limit of 93.5 GeV/c^, no massive particles were found. 
Similar null searches were carried out by Borate et al. [11], the ALEPH experiment 
at LEP. 
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Figure 1: Upper limit of the cross section for pair production of stable charged heavy 
leptons at y/s = 133 to 189 GeV from Acciarri et al. [10]. The hatched area indicates 
the excluded region. The dashed line represents the calculated pair production cross 
section for heavy leptons at y/s = 189 GeV. 
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p. Abrcu et al. [12], the DELPHI experiment at LEP, carried out a more general 
search at 189 GeV for integer charge, stable, massive particles produced in events 
with two or three charged particles. The motivation was to look for the supersym- 
metric partners of the muon and tau, assuming these partners to be stable. No 
signal was found in the search mass range of 2 to 80 GeV/c^, but in this case the 
significance of the search depends upon the assumed cross section. 

Searches for integer charge, stable, massive particles at the Tevatron reach to 
larger masses than searches at LEP2. Abe et al. [13] looked for pair production of 
particles with masses up to the order of 500 GeV/c^. No evidence for such particles 
was found [13]. However, proton- antiproton production cross section calculations 
do not have the certainty of electron-positron electromagnetic pair production cross 
sections, hence the significance of the null results depends upon the model used for 
the production cross section, Fig. 2. 




Figure 2: The cross section upper limits (95% C.L.) for the pair production of 
stable charged fermionic colored particles for Q = 1, Q = | and Q — \ from Abe 
et al. [13]using the Tevatron. Also shown are the theoretical cross sections for the 
production of fermionic color sextets(6), octets(8) and decuplets(lO). 
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2.3 Cosmic rays, halo particles, and searches for integer 
charge, stable, massive particles 

Before discussing experimental limits on X particles in cosmic rays and the galactic 
halo, we turn to some preliminary considerations. First, as noted in Sec. 2.1, much 
of the interest in X particles has been generated in the hope of identifying the X as 
exotic dark matter. This simphfies the search because, as discussed in Sec. 1.4, the 
conditions for the X particles to reside in the halo (a dark matter requirement) are 
that they have a large mass, Mx ~ 10^ GeV/c^, and that their velocity be close to [3, 
4] /3 ~ 10^'^. Therefore one looks for slow, charged particles whose velocities remain 
constant as they pass through the detector. In contrast, a general search for a rare 
X in the cosmic radiation must measure a combination of particle parameters, such 
as mass, velocity, and charge, in order to rule out the large background of naturally 
occurring isotopes. Ultimately, however, it may be impossible to distinguish the 
origin of a detected X particle: it may come from a relic population circulating 
in the halo, a primary cosmic ray made in an astrophysical source or a secondary 
cosmic ray created in a high-energy collision in the earth's atmosphere. Thus, in 
our discussion below we will not distinguish between these but will collectively refer 
to them as cosmic ray particles. 

A second consideration pertains to the identification of X as it passes through 
the detector. Over their history the X"*" and X^ particles may take different guises. 
For example, the X+ could acquire an electron to form an isotope of hydrogen X+e~. 
Similarly, the X~ could have bound to either a proton to form a neutral composite, 
X~p"'", or to an alpha particle which later captures an electron to form a heavy 
hydrogen-like atom, X~a++e~, similar to the X+e~. The X+e~ and X^a++e~ 
are collectively denoted as CHAMPS by De Rujula et al. [3] whereas the X"p+ 
are termed neutraCHAMPS. On entry into the atmosphere, the neutraCHAMPS 
very quickly exchange their proton for ^^N. Thus only satellite or balloon searches 
have to worry about missing a neutraCHAMP because of lack of a signal and, 
conversely, probably only these experiments can distinguish neutraCHAMPS from 
charged CHAMPS by observing an exchange occurring in their detector (e.g.. Bar- 
wick et al. [14], discussed below in Sec. 2.3.1). 

Finally, we consider the upper bound on the mass of X. As discussed in Sec. 1.3, 
the upper bound on the mass of any particle that is produced in the early universe, 
regardless of its interactions, depends on our understanding of the physics of this 
epoch. Although the accepted view holds that Mx < 10^ GeV/c^, a number of recent 
works, for example Chung et al. [15], have described scenarios which allow this limit 
to be extended, even to the GUT scale. Chung et al. [15] resurrect the possibihty 
that not only could dark matter particles be supermassive but that they could also be 
charged. We find, in fact, that the latter possibility can be ruled out by considering 
the results of the MACRO experiment. We discuss this below in Sec. 2.3.2. 
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2.3.1 Balloon and satellite search limits 

Although no satellite or balloon experiment was designed to explicitly search for 
CHAMPS or neutraCHAMPS, limits on their abundances have been obtained us- 
ing data from three such experiments that were originally intended for cosmic ray 
studies. In all three cases the data were obtained and reanalyzed with the intent 
to test the charged dark matter hypothesis proposed by Refs. 3, 4. Prom Eq. 4 
we note that for low mass values, Mx < 10^ GeV/c^, the dark matter flux at the 
top of the atmosphere exceeds that of cosmic rays. This fact is utilized in each of 
these studies to rule out X masses for which the expected dark matter flux is higher 
than observed. Note that the range of excluded X masses also has a lower hmit. 
This is because at small enough X masses the kinetic energy of the particle is no 
longer sufficient for it to penetrate both the solar wind and the minimum number 
of detector elements required by the analysis selection criteria. 

The first limit on the X abundances in cosmic rays [4] was set by using data 
from track-etch detectors employed in the University of Chicago experiment on the 
Pioneer spacecraft. These detectors are ideal for searching for X particles in cosmic 
rays because their restricted sensitivity to particles with charge to velocity ratios 
oi Z/f3 > 14, for (3 > 10~^ naturally eliminates the dominant background of fast 
protons and alpha particles. The response of these detectors to CHAMPS with 
P 10~^ is estimated [4] to be comparable to that of cosmic rays with Z/P 16. 
Ref. 4 also notes however that cahbrations of track-etch detectors are hmited to 
charged particles with f3 > 10~^. Nevertheless, with a number of caveats, dark 
matter being entirely composed of CHAMPS of mass Mx < 10^ GeV/c^ is ruled 
out. However, no limits were obtained on the neutraCHAMP abundance. 

The second analysis, Barwick et al. [14] employed track-etch data from a balloon 
flight and set bounds on both CHAMP and neutraCHAMP masses. This more 
careful analysis includes detailed calculations of CHAMP energy loss, which appear 
to confirm the sensitivity of track-etch detectors to low velocity massive charged 
particles. Barwick et al. [14] rule out dark matter CHAMPS in the mass range 
3.5 X 10^ <Mx < 8.6 X 10^ GeV/c^. In addition they also rule out neutraCHAMPS 
as dark matter for 10^ < Mx < 4 x 10=^ GeV/c^. 

The last analysis, by Snowden-Ifft et al. [16] uses data from the IMP 8 satellite 
to rule out CHAMP dark matter in the mass range 2.4 x 10^ < Mx < 5.6 x 10^ 
GeV/c^. Limits from all three analyses are shown in Fig. 3. 
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Figure 3: Observed upper limits on the flux of CHAMPS (charged massive particles) 
versus the CHAMP mass, Mx- The diagonal line gives the predicted dark matter 
flux as a function of Mx- The excluded, shaded region on the left comes from the 
observations discussed in Sec. 2.3.1. The excluded, shaded regions marked a, b and 
c come from (a) Eq. 7, (b) Eq. 5, and (c) Eq. 6 in Sec. 2.3.2. 



2.3.2 The MACRO experiment secirches for integer charge, stable, mcis- 
sive peirticles 

Of the ground- and underground-based experiments, two place very strong limits 
on the flux of X particles in cosmic rays: MACRO, based deep underground and a 

precursor to MACRO, placed on the surface of the earth. 

The MACRO detector [17], based deep underground in Gran Sasso, Italy, was 
primarily designed to search for magnetic monopoles. In order for the X particle 
to be detected by MACRO, which has a minimum overburden of 3300 ni water 
equivalent, it must be highly penetrating and, therefore, must have a large kinetic 
energy. This places a lower limit on the mass of X. For example, if X has a typical 
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halo particle velocity of /9 ~ 10^"^ at the surface of the earth, it must have a minimum 
mass of Mx,min ~ 10^^ GeV/c^ to reach the MACRO detector. Here we have used 
a mean energy loss [3, 18] of dE/dx ^ 150 MeV/g/cm^ vahd over /3 ~ 10~^ to 
10~^. The MACRO detector consists of a series of subdetectors that are used both 
individually and cumulatively to set flux limits for monopoles and a number of 
other rare exotic particles [19]. Although limits on integer charge particles, such 
as CHAMPS, are not quoted explicitly in Ambrosio et al. [19] they should not 
be too different from their limits on the monopole flux obtained using the liquid 
scintillator subdetector system. Both CHAMPS and monopoles are slow, highly 
ionizing, and penetrating and should give similar responses from scintillators. In 
fact, MACRO was optimized in part by utilizing the measured light yield of slow, 
j3 ~ 2.5 X 10~^ protons in scintillators [20], which should be comparable to that 
of CHAMPS and is estimated to be larger than that due to bare GUT monopoles 
for 10~^ < /9 < 2 X 10~^. The triggers that MACRO employs for monopoles in 
the velocity range 10~^ < /9 < 10~^ should also have comparable efficiency for 
triggering CHAMPS. Therefore, we use the Ambrosio et al. [19] scintillator derived 
90% C.L. upper limit on monopole flux in the velocity range j3 ~ 10^"' to 10"^ as a 
conservative bound on the flux of X particles with mass Mx > 10^^ GeV/c^ in the 
same velocity range: 

$ < 4 X 10-^^cm"^sr"^s"\ (5) 
For faster particles, /3 ~ 1.2 x 10~^ to 10~^, a better limit is given: 

$ < 2.6 X 10"^*^cm-2sr-^s"^ (6) 

which still overlaps with the high velocity tail of the expected dark matter flux. 
In this case, however, the mass reach is lowered to as much as Mx > 10^ GeV/c^ 
because of the higher kinetic energies of the particles. Both of these limits are shown 
in Fig. 3 together with the flux expected from dark matter. 

As mentioned earlier, this important limit on the flux of X particles pertains 
directly to the recent proposal by Chung et al. [15] for supermassive dark matter 
candidates. They alluded to the possibility that such particles could be charged. 
However, if Mx > 10^^ GeV/c^, the above limits surely rule out this possibihty. 

Results from a prcciu^sor experiment to MACRO placed at the surface of earth 
[21], also with the intent to search for magnetic monopoles, can be used to place 
limits on X particles in the mass range Mx > 10^ GeV/c^. Below this mass a 
particle with velocity (3 fa 10~^ at the top of the atmosphere ranges out before 
hitting the ground. Their 90% C.L. upper flux limit for particles with velocity 
5 X 10-^ < /3 < 2.7 X 10-^ is: 

$ < 4.7 X 10-^^cm-V-^s-^ (7) 



13 



As shown in Fig. 3 this hmit, together with those from MACRO and the satelhte 
and balloon limits (Sec. 2.3.1), essentially close the fully allowed mass range for 
charged dark matter of the type considered by Chung et al. [15] 

2.4 Searches in ordinary matter for integer charge, stable, 
massive particles 



2.4.1 Searches for massive particles in water 

It is an attractive idea that a stable, massive, positively charged particle, X"*", falling 
onto and through oceans and lakes, will form the heavy water molecule HXO. Table 1 
lists the four experiments that have searched for HXO. No evidence was found for 
an X"*" in the mass ranges given in the table. The upper limits on the concentration 
of HXO in H2O are given in Fig. 4. 
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Figure 4: Observed upper limits on nHxo/''^H20j the concentration of HXO in H2O 
versus the mass of a unit charge particle X. Limits are from Vekerk et al. [24], 
Hemmick et al. [23] and Smith et al. [22]. 
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Table 1: Searches for stable, massive, positively charged particles in water. The 
upper limits on the concentration of HXO in H2O are given in Fig. 4. 



Concentration Method 


Detection Method 


Experiment 


electrolysis 


mass spectrometer 


Smith et al. [22] 


see paper 


mass spectrometer 


Hemmick et al. [23] 


ultracentifuge 


spectral analysis 


Verkerk et al. [24] 


none 


mass spectrometer 


Yamagata et al. [25] 



2.4.2 Seetrches for massive peirticles in isotopes 

A massive particle with negative integer charge, X~, can electromagnetically bind 
with a very small Bohr orbit to a nucleus [23]. If the nucleus has charge Z, this gives 
an anomalously heavy, usually called superheavy, isotope with charge Z-1. There 
have been several searches for such isotopes. 

Hemmick et al. [23] searched for superheavy isotopes of lithium, beryllium, 
boron, carbon, oxygen and fluorine in a mass range up to about 10^ GeV/c^. No 
superheavy isotopes were found with the upper limits on the number of superheavy 
isotopes per nucleon ranging from 10~^° to 10"^"^. Similar upper limits were found 
by Turkevich et al. [26] for carbon and Dick et al. [27] for sodium. 

2.4.3 Compcirison of bulk matter limits with cosmic ray limits 

The relative fraction of X particles in the halo versus condensed matter in the 
disk (e.g., stars, planets, etc.) is believed to be large if the mass of X is much larger 
than 10^ to 10^ GeV/c^ (Sec. 1.4). This, however, does not mean that a significant 
enrichment could not have occurred in the earth from stopping X particles over many 
years. Therefore, one must weigh the benefits of large throughput experiments in 
bulk matter, which may contain very small abundances of X particles, over cosmic 
ray experiments which have small areas and exposure times but probe the primary 
X particle flux. For this purpose it would be useful to convert the upper bounds 
on the flux of X particles from Sec. 2.3 into an enrichment of ordinary matter on 
earth. Here we describe a simple, illustrative model that does this. 

The assumptions used in the model are extremely simple. The enrichment of 
earth by X particles is assumed to occur continuously over its age, 4 Gyr, during 
which time geophysical or geochemical processes are ignored. The flux of X particles 
and their corresponding masses are taken from the best limits set by experiments 
described in Sec. 2.3 (see Fig. 3). The range of these particles in the earth is 



15 



calculated assuming a constant velocity of f3 ^ 10^''' upon entry, and a mean energy 
loss of dE/dx = 150 MeV/g/cm^. The particles are assumed to stop uniformly 
down to a depth given by their range. Finally, the fraction of X particles in ordinary 
matter, nx/n^, is calculated at the surface and the result is compared in Fig. 5 to 
the results of bulk matter searches given in Fig. 4. 




Figure 5: Upper limits on the fraction of X particles in ordinary matter, hx/un 
where X is a unit charge particle. The left excluded region is from Fig. 4. The right 
excluded region is calculated from the upper limits on the flux of X particles at the 
earth's surface using the model in Sec. 2.4.3. This right exclusion region assumes 
that all X particles in the earth come from the flux of X particles at the earth's 
surface. 



2.5 Limits from astrophysics on the existence of integer 
charge, stable, massive particles 

A number of constraints can be placed on the nature of integer charge massive stable 
particles from astrophysics and cosmology. One cosmological constraint, the upper 
limit on the mass of a particle produced in the early universe, Mx < 10^ GeV/c^, 
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has already been discussed in Sec. 1.3. Another constraint, discussed in Sec. 2.1, is 
that the X particle density cannot exceed the upper limit on the density of exotic 
dark matter which is approximately 30% of the energy density of the universe. The 
local density of this component is believed to be of order 0.3 GeV/c^ which sets 
an upper limit on the flux on X particles at earth (see Eq. 4). Lower bounds on 
Mx have also been set based on theories of galactogcncsis and assuming that X 
comprises 100% of dark matter as given in Ref. 3,4. However, as discussed in Sec. 
1.4, their limits, 2 x 10^ GeV/c^ and 10^ GeV/c^, disagree by a large margin. If 
we adopt the conservatively calculated higher limit of Dimopoulos et al. [4], and 
believe the upper bound Mx < 10® GeV/c^ set from cosmology, this essentially 
rules out the CHAMP dark matter scenario. Below we outline two more examples 
of astrophysical bounds to give a flavor of this type of inference. 

Ghivukula et al. [28] rule out GHAMP dark matter based on measurements of 
cooling rates of clouds of hydrogen gas in the interstellar medium. Their argu- 
ment relies on noting that in equilibrium the heating of the interstellar gas from 
interactions with CHAMPS cannot exceed the cooling rates. This places an upper 
bound on the strength of interactions between halo dark matter and the interstellar 
medium. By estimating the elastic GHAMP-hydrogen interaction cross section and 
using the local dark matter density, they find that Mx > 10^ GeV/c^. This limit 
only applies to CHAMPS (i.e., X+e" and X-a++c") and not to neutraCHAMPS. 

Gould et al. [29] also rule out CHAMP dark matter but in a large mass range, 
10^ < Mx < 10^® GeV/c^. They argue that the hypothetical dark matter CHAMPS 
(X+e~ and X~Q;++e~) would accumulate in protostellar clouds with a sufficient 
concentration to disturb latter phases of stellar evolution. Specifically, in the case 
where a neutron star is the endpoint of a star's life, the CHAMPS drift to its center 
and form a black hole which quickly destroys its host via accretion. For dark matter 
CHAMPS in the mass range given above, Gould et al. calculate that the destruction 
of the neutron star would occur on time scales less than 10 yr, which is much shorter 
than the ages of old neutron stars (>1 Gyr). 

3 Searches for Neutral, Stable, Massive Particles 

Weakly interacting massive particles (WIMPS) are considered one of the leading can- 
didates for cold dark matter [6,30,31]. Of the various types of WIMPS, the lightest 
supersymmetric particle (LSP, e.g., neutralino) with mass between 1 GeV/c^ and 1 
TeV/c^ and interaction cross section with ordinary matter a < a^eak, is amongst the 
favored [31]. In Sec. 3.1 we briefly review the experimental status of dark matter 
WIMP searches and show the current best limits. In Sec. 3.2 we discuss another 
candidate for dark matter, strongly interacting massive particles (SIMPS), which, 
except for a small region of mass versus cross section space, have been ruled out. 
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3.1 Searches for WIMPS - weakly interacting massive par- 
ticles 



There have been numerous reviews on the subject of WIMPS [6,30-35]and descrip- 
tions of various experiments and their status [32,36]. Here we discuss some general 
features of dark matter WIMP searches and briefly review the most sensitive cur- 
rent experimental limits. Ellis et al. [37] have reviewed accelerator searches and 
limits. There are two basic methods for detecting dark matter WIMPS: the direct 
and the indirect search methods. In the direct method (Sec. 3.1.1) the signal of 
a WIMP interacting in the detector (e.g., from WIMP-nucleon elastic collision) is 
directly measured. In the indirect method (Sec. 3.1.2) it is the products of WIMP- 
antiWIMP annihilation (e.g., gamma rays, positrons, anti-protons, or neutrinos) 
which are detected. Presumed sources can be in the galactic halo or in specific sites 
where WIMPS can gravitationally accumulate over time, such as at the center of 
the Galaxy, earth or sun. 

3.1.1 Direct searches for dark matter WIMPS 

In this category there are currently over twenty experiments that are in various 
stages of planning, construction, or running [31,38]. One experiment, DAMA, has 
in fact claimed a positive signal, which they find to be consistent with a WIMP 
mass of 52l8° GeV/c^ at 4o" CL. However, a second experiment, the Cryogenic 
Dark Matter Search (CDMS), finds a null result even though their sensitivity is 
close to that of DAMA's. Below we review the results of these experiments but, 
first, we discuss some general features of direct dark matter WIMP searches. 

The differential interaction rate of WIMPS in a detector, in units of rate per 
unit energy per unit detector mass, is given by [30]: 



Here pdm — 0.3 GeV/cm^ is the local density of dark matter (the assumed WIMP 
density), mw is the WIMP mass, cr is the WIMP-nucleon elastic interaction cross 
section, ^ is the WIMP-nucleon reduced mass, F{Q) is a nuclear form factor, Vmin 
is the minimum detectable WIMP velocity, and f{v) is the distribution of WIMP 
velocities relative to the detector. Also 



where Ethr is the detector energy threshold and is the target nucleus mass. Note 
that the form factor F{Q) depends on WIMP-nucleus coupling, which can be spin- 
dependent or spin-independent. In general the spin-independent coupling dominates 






(9) 
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and it is from results of experiments probing this interaction that the most sensitive 
bounds exist. See Fig. 5 of Ref. 34 for hmits from recent spin-dependent results. 

The unknowns in Eq. 8 which are sought are rriw and a. The astrophysical 
quantities pou and f{v) are determined from models of the galactic halo in which 
the WIMPS are assumed to be distributed in a non-rotating isothermal sphere with 
a velocity distribution close to Maxwcllian [6,30] (see Sec. 2.1). However, see Evans 
et al. [39], Green [40] and UUio and Kamionkowski [41] for discussions of deviations 
from the standard halo model. In this scenario the sun moves through the WIMP 
"wind" with an average velocity [30] of 220 km/s, whilst the earth's velocity relative 
to the wind varies by a few percent as it revolves around the sun. This variation 
should result in a seasonal (sinusoidal) modulation of the event rate together with 
a variation of the energy spectrum. The event rate should be a maximum when 
the earth's motion is aligned with the sun's (June) and a minimum when it is anti- 
aligned (December). The effect, however, is only of order several percent and to be 
seen it requires a large detector to acquire the necessary statistics. Nevertheless, if 
seen, it is one of several important criteria that can be used to separate signal from 
background [6]. These seasonal effects on the rate and spectrum are, in fact, what 
DAMA observes and are the main basis for their claim. 

To understand detector requirements it is important to know the rates expected 
from dark matter WIMPS. Typical rates predicted from supersymmctry models 
range from 10"^ to 10 events/kg/day, for a deposition of energy between 1 and 100 
keV in the detector [31]. If one includes cosmological constraints, such as requiring 
WIMPS to fill the dark matter budget, then one gets rates < 10~^ events/kg/day, 
which are about an order of magnitude below current detector sensitivities. Given 
the rarity of WIMP interactions it is important to understand and use specific signa- 
tures that can be used to eliminate backgrounds and bolster claims [6,30,31]. Two 
of these, the modulations of event rates and corresponding variations of the energy 
spectrum, have already been mentioned. Another important WIMP signature can 
be obtained if the directional information of the nuclear recoil in the detector is 
available. The nuclear recoil is anisotropic due to the seasonal modulation of the 
flux discussed above. However, the recoil signature is larger and requires fewer 
statistics to resolve [42]. The DRIFT experiment is the only one currently running 
which employs this signature [43]. Also important are signatures which specifically 
distinguish the WIMP signal from cosmic ray backgrounds and radioactivity back- 
grounds which are by far the largest and most insidious. The reduction of this 
background requires using a combination of radio-pure detector elements, shielding, 
and the ability to distinguish WIMP interactions in the detector from those of the 
background. The WIMPS (and neutrons) deposit a much larger fraction of their 
energy as heat (via nuclear recoils) whereas interactions from backgrounds such as 
gammas, electrons, and alpha particles result in ionization from electron recoils. The 
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CDMS experiment, for example, measures both heat and ionization deposition in 
their detectors, using the ratio of the two to segregate the background events from 
the neutron or possible WIMP events. Other techniques for background reduction 
are given by [32,38,41,44]. Other signatures of WIMP detection can be found in 
Ref. 6. 

The current experimental state is dominated by the recent discovery claim from 
the DAMA experiment [45] at the Gran Sasso National Laboratory of I.N.F.N. 
DAMA employs large radiopure Nal crystal scintillators as their detectors. They 
use 9 Nal crystals, 9.7 kg each, allowing them to obtain 57,986 kg-day worth of 
statistics spanning four different yearly cycles (which are also reported individually 
as DAMA/NaI-1 through DAMA/NaI-4). From calibration data they determined 
that nuclear recoils from neutrons produced faster signals in the Nal detectors than 
electron recoil pulses from gammas or electrons. This timing difference is used in 
a pulse shape analysis for background discrimination after which their background 
rate is lowered to 0.5-1 events/keV/kg/day in the 2-3 keV part of the spectrum, and 
about 2 events/keV/kg/day in the 3-20 kcV region [32]. Their claim of finding a 
WIMP {mw = 52l^° GeV/c^ and = 7.2;[J:^ x 10"^ pb, at Aa C.L.) is based on 
their observing a modulation in their background-subtracted rates consistent with 
that due to the seasonal modulation signature discussed above. Figure 2 of Ref. 45 
shows the modulation signature. Here f = -^^is the fraction of the dark matter 

^ ^ PDM 

density made up of WIMPS; it appears often in the dark matter literature. Thus 
the cross section is uncertain by the factor ^. 

This result, however, has met with some skepticism and criticism, most of which 
is centered on the shape of the background-subtracted spectrum [31,46] . In addition, 
raising further questions [31, 77], other Nal experiments employing similar pulse 
shape discrimination have encountered a class of unexplained events with pulse times 
close to, or even shorter than, those expected from neutrons or WIMPS. Finally, the 
CDMS experiment, with sensitivity similar to that of DAMA, has recently released 
results that are inconsistent with DAMA's allowed ruw — ^cr region [38,44]. 

The CDMS experiment [38,44] is currently housed in the Stanford Underground 
Facility. CDMS has the ability to measure both phonons and ionization on an event- 
by-event basis. Using the ratio of energy-loss in phonons to ionization, which for 
nuclear recoils is lower by a factor of 3 than for electron recoils, the CDMS exper- 
imenters are able to achieve a background reduction of greater than 99%. This, 
together with the use of radiopure detector elements and care in shielding, has 
allowed CDMS to reach background rates below 0.1 events/keV/kg/day in the 10- 
20 keV region of the spectrum [47]. Thus, although the small size of the CDMS 
apparatus limits exposure to only 10.6 kg-days, the ability to efficiently reject back- 
grounds makes the CDMS experiment competitive. Prom the data set achieved in 
this exposure, they report 13 events resulting from nuclear recoils between 10 and 
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100 keV [38]. Although this rate is similar to that expected from DAMA's claimed 
signal region, the CDMS analysis of these events finds the events to be consistent 
with neutrons [38] . Thus the findings of the CDMS experimenters rule out DAMA's 
signal region at greater than 75% CL. In addition, the CDMS experimenters find 
their data to be incompatible with the DAMA modulation signal shown in Fig. 2 of 
Ref. 45 at greater than 99.98%. The CDMS and DAMA results are shown in Fig. 
6. 

The peculiar situation existing today in the field of dark matter searches surely 
will not last. CDMS will soon move to the Soudan Mine where the background rate 
should go down to 10~^ to 10~^ events/keV/kg/day [47]. The resulting increase in 
the expected sensitivity of CDMS at Soudan, shown in Fig. 6, will allow it to probe 
the mw—^cr space some 2 orders of magnitude below DAMA's allowed region. Other 
experiments coming on line will have similar or even better reach and together they 
will soon be in a position to verify or rule out the DAMA result. 




Figure 6: Current observations and limits on the WIMP-nucleon spin-independent 
cross sections as a function of mass. See the text for explanation of the curves and 
areas. 
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In Fig. 6, the top three curves are upper hmits from: (sohd) the Heidelberg- 
Moscow experiment [48], (dash-dash) the 1996 DAMA experiment [45] and (dot- 
dash) the CDMS experiment [38, 44] . The hmit from the Heidelberg- Moscow ex- 
periment [48] is shown because it has achieved the lowest background rates directly 
from raw data. The dark shaded area is the observed 3a signal from the DAMA 
experiment. The light shaded area shows the theoretical parameter space for various 
supersymmetric models for WIMPS. The dot-dot curve shows the projected upper 
limit that may be reached by the CDMS experiment when it is in the Soudan mine. 
The figure was obtained from the limit plots web tool provided by R. Gaiskell and 
V. Mandic at http://cdms.berkeley.edu/limitplots. At this web site, combinations 
of current and projected future experimental limits for both spin-dependent and 
spin-independent cross sections can be viewed. 

3.1.2 Indirect secirches for deirk matter WIMPS 

The detection of dark matter through its annihilation signal from various astrophys- 
ical sites is the concept behind a number of experimental programs [33,47]. WIMPS 
should accumulate gravitationally in the center of the sun, and the earth, increas- 
ing the annihilation rate relative to that from WIMPS in the halo. WIMPS could 
also accumulate in the center of the galaxy [49]. In all three cases the resulting 
enhancement in the directional anisotropy from the decay products can be utilized 
as a WIMP signature. WIMP annihilation in the halo, although carrying much less 
directional information, may produce excesses of gamma rays, positrons or anti- 
protons above that expected from standard astrophysical processes [50-53]. There 
have been a number of reports of features in positron [53] or anti-proton spectra in 
cosmic rays [54] which have been interpreted as coming from WIMP annihilation in 
the halo, but none has been confirmed [33]. 

Gamma rays and neutrinos provide a somewhat clearer signal because, unlike 
cosmic rays, they carry directional information. Furthermore, for gamma rays the 
spectrum can be probed for features from monoenergetic photons produced in WIMP 
annihilation in various scenarios. For example, if WW 77 then ?a mw] 
see Ref. 50 for further examples. In fact, simulations shown in Ref. 33 of the 
sensitivities of the next generation ground and space-based gamma ray detectors to 
WIMP annihilation products 77 and Z7 indicate that competitive constraints could 
be placed in supersymmetric model parameter space. 

The possible accumulation and subsequent trapping of WIMPS in the sun and 
earth are the best hope for using neutrinos as probes of WIMP annihilation [33] . 
High energy neutrinos, the final products in a majority of WIMP annihilation 
schemes, can provide a clear WIMP signal if observed coming from the sun and 
the earth. Calculations [33,47] of event rates in detectors from these energetic neu- 
trinos indicate that a detector of area 1 km^ could begin probing interesting regions 
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of supersymmetric parameter space. A number of large neutrino detectors (e.g., 
AMANDA, Baksan, MACRO, and Super-Kamiokande) have reported flux limits 
from point sources, although none has shown an excess from any of those consid- 
ered above. Reference 55 reviews recent results. The advent of next generation large 
neutrino detectors should provide competitive and complementary limits to direct 
searches for WIMP dark matter. 

3.2 SIMPS 

In the mid-1980's it was noted that while attention was centered on WIMPS, strongly 
interacting massive particles (SIMPS) could also satisfy the dark matter 
hypothesis [56,57]. In fact, their interaction cross sections with matter he anywhere 
from a > (Tweak up to a barn. Numerous candidates were suggested [56,57], including 
both elementary particles and bound states, and it was suggested that they could 
be detected using detectors similar to those used to search for WIMPS [56]. Soon 
after the CHAMPS hypothesis (see Sec. 2), Starkman et al. [58] applying a similar 
analysis, calculated existing constraints on the cross sections of dark matter SIMPS 
for masses greater than about 1 GeV/c^. They found a few interesting open mass 
ranges: m < 10^ GeV/c^ with cross sections greater than 10^^^ cm^; 10^ < m < 10^ 
GeV/c^ with cross sections greater than 10^^^ cm^; and m > 10^^ GeV/c^ with a 
very broad range of cross sections that extends to below 10^^^ cm^. In addition, for a 
number of these regions, the viability of SIMPS as dark matter depends on whether 
or not there is SIMP-antiSIMP symmetry in nature. See Figs. 1-3 of Starkman et 
al. [58] for details. Below we briefly summarize the experimental efforts that have 
been undertaken to close the SIMP dark matter windows. 

A combination of balloon-borne, ground-based, and underground experiments 
have essentially closed the 10^ < m < 10^ GeV/c^ mass range for dark matter 
SIMPS [59-61]. In addition, the underground experiment (DAMA), which is ideal 
for searching for rare penetrating particles, has closed much of the high mass m > 
10^° GeV/c^, a > 10^^° cm^, region [61]. Mohapatra and Nusinov [62] show that 
other WIMP searches, together with heavy isotope searches, also rule out the mass 
range from a few GeV/c^ to a few TeV/c^ for dark matter SIMPS. 

Although the above constraints have effectively closed off much of the interesting 
phase space for dark matter SIMPS, interest in them has not died out. Teplitz et 
al. [63] have recently proposed laboratory experiments for searching for anomalous 
heavy nuclei such as gold. Here, however, they specifically relax the constraint 
that SIMPS satisfy dark matter abundances. Also, SIMP motivations for dark 
matter have been revived in the form of strongly self-interacting particles which have 
been suggested as a means to circumvent astrophysical and cosmological problems 
perceived with the cold dark matter model [64]. 
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4 Searches for Fractional Charge, Stable Particles 



In this section we ignore the restriction to massive particles. The discovery of a 

stable, fractional charge particle of any mass would be of great importance. 

In the 1970's and 1980's there were four reviews of searches for fractional charge 
particles. General reviews, including accelerator and cosmic ray searches, were given 
by Jones [65] and Lyons [66]. Smith [67] primarily reviewed searches in terrestrial 
materials. MarineUi and Morpurgo [68] reviewed the magnetic levitation method. 
These reviews provide much information on search technologies. Our review is 
devoted to the results of the searches of the 1990's with older results included when 
still relevant. 

We give all charges in units of the magnitude of the electron charge, e. 

4.1 Accelerator searches for fractional charge, stable parti- 
cles 

High energy, high intensity, electron-positron and proton-antiproton colliders are 
ideal instruments for searches for fractional charge particles. Surprisingly with re- 
spect to published results, full use has not been made of these opportunities. One 
limitation is that the experiments, or at least the published analyses, have con- 
centrated on q = |e, |e, and |e. However in some cases the search may be more 
general. 

4.1.1 Seairches using electron-positron colliders 

The study of Z° decays provides a direct and an indirect way to search for fractional 
charge particle pair production at the Z^: 

e+ + e- ^ Z° ^ + X"^. (10) 

However, there appear to be no published reports of direct searches at the ^ for 
fractional charge particles produced according to Eq. 10. 

The indirect search method is based upon the concept that if the X pair is 
not detected, then its decay width contributes to the invisible width of the Z°. X 
particles with g < 1 might not be detected because of their smaller than normal 
gas ionization in the track detectors. In units of neutrino width contribution, the 
direct measurement of the invisible width gives the quantity 3.07±0.12. Therefore, 
if fractional charge particles are pair produced but not detected at the Z°, their 
weak interaction coupling must be smaller than that of the neutrino by a factor of 
order 2 or 3. This condition on the weak interaction coupling constant is thus a 
hmit on the existence of a fractional charge particle with mass < 45.5 GeV/c^. 
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Two searches at the Z° using inclusive production have been reported. Buskuhc 
et al. [69], the ALEPH collaboration, searched for particles with q = ^e, |e, or |e 
in the mass range of 8 to 45 GeV/c^. Their 90 % C. L. upper hmits are given 
in terms oi R = o^xx/o'^i^i where o"^^ is the electroweak pair cross section. The 
upper limits on R are in the range of 1 to 3 x 10~^. Akers et al. [70], the OPAL 
collaboration, reported similar null results. 

AckerstafT et al. [71], the Opal collaboration, carried out a LEP search at higher 
energy, 130 to 183 GeV for pair produced X particles with |e. The X mass search 
range was 45 to about 90 GeV/c^. The 95 % C. L. upper hmits on the cross section 
were in the range of 0.05 to 0.2 pb. 

There is an experimental uncertainty when searching at a collider for fractional 
charge particles that have quark-like charges. Is one actually searching for free 
quarks and, if so, what is the interaction cross section of the free quark? If it 
is a strongly interacting cross section then the particle may turn into conventional 
hadrons before its charge can be measured. Guryn et al. [72] considered this problem 
in their 1984 search, e+ + e~ — X+'^ + X~'^, at 29 GeV. They used a detector with 
a thickness of less than 1 % of a conventional hadronic interaction length. No 
fractional charge particles were found , the search range extending to 13 GeV/c^. 

4.1.2 Seairches using proton-antiproton colliders 

The highest energy reported search, p + p — > X+^ + X~^, for fractional charge 
particles is that of Abe et al. [13] carried out at the 1.8 TeV Tevatron. Figure 2 
shows the 95 % G. L. upper limits on the cross section for q — and |e. 

4.1.3 Searches using heavy ion collisions 

In the 1980's there were two accelerator searches using heavy ion collisions based on 
the concept that "processes that involve heavy nuclei might yield fractional charge by 
mechanisms not possible in elementary particle interactions" to quote Lindgren et al. 
[73]. These experimenters examined the product of the collision of 1.9 GeV/nucleon 
^^Fe nuclei with Pb nuclei and gave an upper limit of less than 10~^ fractional charge 
particles produced per Fe-Pb collision. 

In another experiment, Barwick et al. [74] examined the 14 < Z < 18 projectile 
fragments of 1.8 GeV/nucleon "^"Ar collisions searching for fractional charge particles 
with null results. 
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4.2 Cosmic rays, halo particles, and searches for fractional 
charge, stable particles 

The most comprehensive search to date for fractional charge particles impinging on 
the earth has been carried out by Ambrosio et al. [75] using the MACRO experi- 
mental apparatus. They searched for lightly ionizing particles with .25 < /5 < 1.0 
and were sensitive to charges as small as No fractional charge particles were 
found, the 90% C.L. upper limit on the flux is of the order of 

< 2 X 10"^^cm~V"^s"\ (11) 

The limit as a function of q is given in Fig. 7. Recall that the MACRO ex- 
periment lies beneath a mountain with a minimum overburden of 3300 m water 
equivalent. This leads to two limitations on the significance of these flux limits: 
strongly interacting fractional charge particles would not reach MACRO and non- 
strongly interacting particles must have sufficient energy to overcome the dE/dx 
loss. 

Two other searches for fractional charge particles in cosmic rays and halo par- 
ticles were carried out in the 1990's; both report upper limits on the flux only for 
particles with charges of |e and |e. 

Aglietta et al. [76] used the LSD detector in the Mont Blanc tunnel with an 
overburden of 5000 m of water equivalent to find the 90% C.L. upper limits on the 
flux: 

1 2 

0(-e) < 2.3 X 10-^^cm-2sr-V\ 0(-e) < 2.7 x lO-^^cm-^sr-^s'^ (12) 
3 3 

Mori et al. [77] used the Kamiokande II detector with an overburden of 2700 m 
of water equivalent to flnd the 90% C.L. upper limits on the flux: 

0(-e) < 2.1 X lO-^^cm-^sr-VS (t){\e) < 2.3 x lO'^^cm'^sr-^s'^ (13) 
3 3 

Both sets of limits are shown in Fig. 7. 

4.3 Searches in ordinary matter for fractional charge, stable 
particles 

At present there are two methods for searching for fractional charge, stable particles 
in ordinary matter- the levitometer method and the Millikan drop method. In the 
levitometer method, a small object consisting of ordinary matter is magnetically 
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Figure 7: The curve is the observed upper hmit on the flux of fractionally charged 
particles as a function of e/q from the MACRO detector, Ambrosio et al. [75]. 
The range of (3 is 0.25 to 1.0. The points are upper limits from the Kamiokande 
II detector, Mori et al. [77], and the LSD detector, Aglietta et al. [76]. q is the 
fractional charge in units of the electron charge e. Thus q =1/2 correspond to 
e/q =2. 



suspended in vacuum using either ferromagnetism [67, 68] or superconductivity [78] . 
The object often, but not necessarily, is a sphere. The object's mass is in the range 
of 0.03 to 0.1 mg. An electric field is used to drive the object into forced oscillation 
and the charge on the object is measured by the resonant frequency of the oscillation. 

The Millikan drop method goes back to his measurements of the electron charge 
[79]. In the modern development of the method [80-83], liquid drops with diameters 
in the range of 5 to 30 /im fall though air in the presence of an oscillating electric 
field. The drop charge is determined by measuring the drop's terminal velocity 
caused by the electric field. 

Table 2 lists recent sensitive searches for fractional charge particles in ordinary 
matter. All experimenters reported null results except La Rue et al. [78] who re- 
ported finding |e and |e charges in niobium. This 1981 paper produced considerable 
interest in the possible existence of free quarks, but Smith et al. [84] studied about 
four times as much niobium and found no evidence for any fractional charge particles 
in niobium. The La Rue et al. [78] report is not accepted today. 
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Table 2: Searches for fractional charge particles in ordinary matter. All experi- 
menters reported null results except La Rue et al. [78]. See text. There are 6.4x10^° 
nucleons in a milligram. 



Method 


Experiment 


Material 


Sample Mass(mg) 


superconducting levitometer 


LaRue et al. [78] 


niobium 


1.1 


ferromagnetic levitometer 


Marinelli et al. [68] 


iron 


3.7 


ferromagnetic levitometer 


Smith et al. [84] 


niobium 


4.9 


ferromagnetic levitometer 


Jones et al. [86] 


meteorite 


2.8 


liquid drop 


Joyce et al. [80] 


sea water 


.05 


liquid drop 


Savage et al. [87] 


mercury 


2.0 


liquid drop 


Halyo et al. [88] 


silicone oil 


17.4 



The sensitivity of the searches in Table 2 may be estimated by noting that there 
are 6.4 x 10^° nucleons in a miUigram. The largest individual search sample is 17.1 mg 
of silicone oil, Halyo et al. [88] . They report that the concentration of particles with 
fractional charge more than 0.16e from the nearest integer is less than 4.7x10"^^ 
particles per nucleon with 95% confidence. 

As shown in Table 2, most materials used for fractional charge searches have been 
chosen for ease of use: niobium for superconducting levitation, iron for ferromagnetic 
levitation, mercury and oil for the Millikan drop method. These are probably not 
the best choices because fractional charge particles might easily be lost in electrically 
conducting materials or in refined materials such as oils. Lackner and Zweig [89] have 
discussed the chemistry of atoms containing fractional charge and have discussed 
the most suitable materials for fractional charge searches [90]. As summarized by 
Perl and Lee [82] the most suitable materials appear to be meteoritic material from 
asteroids, terrestrial minerals that concentrate rare impurities and perhaps material 
from the moon's surface. 

4.4 Searches for millicharge particles 

All the fractional charge searches so far discussed are limited by their inability to 
detect very small charges, less than about e/10. This is because in all the methods 
the 3(7 or 4(j limits on charge measurement precision are about e/10. A direct search 
for particles with charges in the range of 10~-'^e to 10~^e was carried out by Prinz 
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et al. [5] using a high energy electron beam to try to produce pairs of milhcharge 
particles. The search, Fig. 8 from Prinz et al. [5], ruled out a mass region from 0.1 
MeV/c^ to 100 MeV/c^, the 95% C. L. upper limit on the charge being about 10~^e 
at the low mass boundary and 6xl0~^e at the high mass boundary. Figure 8 also 
shows other limits on the existence of millicharge particles [91-93]. 
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Figure 8: Excluded sections of charge- mass parameter space [5]. The dark central 
region is the area excluded by Prinz et al. [5]. The lighter shaded regions give limits 
derived from other experiments. Other regions are excluded by astrophysical or 
cosmological arguments. The references are: (a) Davidson et al. [91], (b) Golwich 
and Robinett [92], and (c) Davidson and Peskin [93]. 



5 Summary 

There is no confirmed experimental or observational evidence for the existence of 
massive, stable, elementary particles other than the electron and proton. In par- 
ticular, in the special case of dark matter there is no confirmed evidence for dark 
matter being composed of massive particles. 

There are many areas where continued or new searches for massive particles are 
warranted: 
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• With respect to searches using proton-antiproton colhders, it seems to us that 
searches at the Tevatron are incomplete. Of course, thorough searches should 
be carried out when the Large Hadron Collider is operating. 

• With respect to searches using electron-positron colliders, it seems to us that 
there is uncertainty in the completeness of searches at the Z° peak and below 
the Z° peak. Of course, thorough searches should be carried out when an 
Electron- Positron Linear Collider is operating. 

• The various searches for fluxes of massive particles impinging on the earth, be- 
gun as dark matter explorations, should be continued even when the observed 
upper limit on the flux is too small to explain the expected dark matter den- 
sity. There arc two reasons. First, dark matter might consist of several types 
of massive particles. Second, the existence of a stable massive particle is of 
great importance even if it has nothing to do with dark matter. 

• Existing searches in bulk matter for massive particles with fractional or integer 
electric charge have probed to sensitivities in the range of 10^^^ to 10^^^ mas- 
sive particles per nucleon. Improved search technology would allow substantial 
improvement of these sensitivities as well as the extension of the searches to 
matter from meteorites and the moon. 
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